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NOMENCLATURE 


A area  of  melting  surface 

on  defined  by  eq  4 

o0  defined  by  eq  6 

B a constant 

c specific  heat  of  water 

9 gravitational  acceleration 

Gr  Grashof  number,  \gL * |0,.(  fp- f.) 1 1 v1 
/>  average  heat  transfer  coefficient 

A thermal  conductivity  of  water 

Lf  height  of  the  ice  sample 

L latent  heat  of  fusion  of  ice 

M melting  rate 

NuL  Nusselt  number  (/>{.)/*  with  melting 

Nu0  Nusselt  number  without  melting 

Pr  Prandtl  number  a $cn/k 

qf  heat  used  for  melting 

qs  sensible  heat 

R radius 

ReL  Reynolds  number,  ( vpL)lu 

Tt  initial  ice  temperature 

Tm  bulk  water  temperature 

Tm  melting  temperature 

fp  plate  temperature 

ST  temperature  difference  (Tm-Tm) 
v water  velocity 

yx  velocity  along  the  plate 

W amount  of  ice  melted 

x coordinate  along  the  plate 

y coor&mete  perpendicular  to  the  plate 

a thermal  diffusivity  of  water 

0 defined  as  c( Tm~Tm)ILf 

0„  defined  by  eq  1 5 

8 time 

(t  viscosity  of  water 

v kinematic  viscosity 

p density  of  water 

Pi  density  of  ice 

r0  shear  stress  at  the  interface  without  melting 

rw  shear  stress  at  the  interface  with  melting 

X a dummy  variable 
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SUMMARY 

An  experimental  device  was  developed  to  simulate  the  heat  transfer  behavior  of  a rather  large  proto- 
type ice-water  heat  sink.  Since  both  the  laboratory  setup  and  the  prototype  heat  sink  used  the  same 
fluid  medium  (i.e.  ice-water),  it  would  have  been  nearly  impossible  to  provide  the  dynamic  similarity 
needed  to  ensure  that  the  laboratory  work  would  reproduce  heat  transfer  characteristics  of  the  proto- 
type. This  difficulty  was  circumvented  by  installing  belts  on  two  opposite  sides  of  the  square  test  col- 
umn. The  belts  were  designed  and  adjusted  to  move  at  the  same  direction  and  velocity  as  the  water 
flowing  through  the  test  chamber.  In  doing  this,  the  laboratory  setup  could  be  considered  as  a st< 
mental  ring  cut  from  the  proposed  prototype. 

The  experiments  covered  bulk  water  temperatures  ranging  from  1 .11°  to  7.50°C  and  water  velocities 
from  1 .7  to  9.8  mm/s.  The  experimental  results  can  be  expressed  in  terms  of  Nusselt,  Prandtl  and 
Reynolds  numbers  as  NuL/Pr  = 3.275(ReL)°  270  with  a moderate  correlation  coefficient  of  0.843.  The 
results  were  also  compared  with  those  computed  from  the  case  of  laminar  heat  transfer  over  a horizontal 
melting  plate.  It  was  found  that  the  results  were  four  to  five  times  higher  than  those  for  the  horizontal 
plate. 

From  this  study  it  was  also  concluded  that,  for  the  case  of  forced  convective  melting  heat  transfer, 
the  effect  of  buoyancy  force  created  by  maximum  density  at  4°C  on  heat  transfer  appeared  insignifi- 
cant even  at  the  lowest  velocity  (i.e.  1.7  mm/s)  conducted  in  this  investigation.  The  heat  transfer  was 
drastically  reduced  in  the  vicinity  of  about  5.6°C  in  the  case  of  free  convective  melting  heat  transter  for 
a wide  variety  of  melting  geometries. 
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HEAT  TRANSFER  OVER  A VERTICAL 
MELTING  PLATE 

Yin-Ouo  Yen  and  Mark  M.  Hart 


INTRODUCTION 

Problem 

Ice  has  been  considered  as  a heat  dissipating  medi- 
um for  enclosed  generating  systems,  and  considerable 
interest  has  been  devoted  to  investigating  the  feasibil- 
ity of  using  ice  in  heat  sinks  for  underground  hardened 
military  installations.  These  sinks  would  be  used  in 
an  emergency  situation  to  absorb  the  waste  heat  cre- 
ated by  an  underground  power  generating  station  and 
its  associated  support  facilities.  During  most  of  a 
power  plant’s  operation,  heat  rejection  is  accomplished 
at  the  ground  surface  with  typical  cooling  towers  and 
water  reservoirs,  but  a system  that  is  capable  of  ab- 
sorbing all  the  waste  heat  for  a short  period  of  time 
must  be  available  at  a moment's  notice.  A typical  in- 
stallation would  consist  of  a nuclear  steam  turbine 
generator,  a steam  condenser  and  three  ice-water  heat 
sinks.  One  concept  in  particular  uses  cylindrical  ice 
columns  of  33.S  m height  and  19.8  m diameter,  with 
a circulating  coolant  rate  of  9.08x  1 0s  kg/h  and  a 
minimum  anticipated  heat  rejection  rate  of  6.73x  10* 
kcal/h  (see  reference  2). 

Brown  and  Quinn2  studied  such  a heat  sink  by  geo- 
metrically scaling  their  experimental  setup  to  the  actual 
33.5-x  19.8-m  prototype  with  a cylindrical  block  of  ice 
1 .83  m in  height  and  1 .22  m in  diameter.  Since  the 
same  fluid  (ice-water)  was  used  for  both  the  model  and 
the  prototype  sink,  it  was  impossible  to  devise  a model 
with  total  similarity  to  the  prototype.  They  developed 
a scale  model  taking  into  consideration  temperature, 
heat  transfer  and  time  factors,  in  addition  to  geo- 
metrical and  performance  characteristics,  as  a basis 
of  analysis  for  the  prototype.  The  results  of  the 
experimental  model  studies  were  compared  with  those 
predicted  from  numerical  computation  of  the  system. 
Favorable  comparison  was  obtained,  providing  confi- 
dence in  using  model  study  results  to  calculate  the 
prototype  performance.  Stubstad  and  Quinn* 


investigated  this  problem  further  by  conducting  a 
larger  scale  experiment  using  an  ice  cylinder  3.0  m 
in  height  and  1 .83  m in  diameter.  Both  studies  have 
provided  a quantitative  description  of  the  melting 
pattern  as  a function  of  operating  time. 

In  a more  recent  study,  Grande*  presented  an  analy- 
sis and  development  of  the  conceptual  design  of  an  ice- 
water  heat  sink  system  to  accommodate  the  closed  cycle 
operation  of  a 1 500-kW  nuclear  power  plant  in  a hard- 
ened underground  installation.  This  study  provided  a 
critical  review  of  previous  work  on  ice-water  heal  sinks. 
It  also  extended  the  analytical  description  of  the  sink's 
thermal  performance  and  assessed  the  available  options 
for  varying  the  flow  path  of  coolant  water  within  the 
sink.  Computer  programs  were  developed  for  predict- 
ing the  time  histories  of  heat  sinks  in  conjunction  with 
power  plants,  and  a practical  sink  design  was  developed 
based  on  the  power  plant  performance  constraints. 

The  present  study  was  conducted  in  order  to  simu- 
late the  exact  flow  conditions  found  within  the  proto- 
type. An  inherent  problem  with  scaling  a flow  situation 
such  as  this  is  that  the  flow  dynamics  of  a scale  model 
cannot  be  exactly  reproduced  in  the  prototype.  To  cir- 
cumvent this  difficulty,  the  present  experimental  setup 
is  designed  in  such  a manner  that  it  is  essentially  a seg- 
ment of  the  prototype.  A melting  ice  plate  is  used  to 
simulate  the  vertical  ice  face  of  the  column  while  the 
back  wall  of  the  test  tank  acts  as  the  containment  wall 
of  the  prototype.  Since  a segment  within  the  annular 
space  of  a prototype  would  be  bounded  on  both  sides 
by  water  flowing  downward  at  the  same  rate,  the  de- 
sign of  the  experimental  setup  includes  continuous  belts 
at  the  opposing  sides  of  the  test  tank  moving  at  the 
same  velocity  as  the  water. 

General  review  of  the  subject 

Most  of  the  previous  work  involving  melting  of  ice 
was  on  free  convective  heat  transfer.  This  work  is  com- 
plicated by  water’s  intriguing  feature  of  reaching  a 
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maximum  density  at  about  4°C.  Since  existence  of 
maximum  density  implies  that  the  thermal  expansion 
coefficient  0 changes  its  sign,  the  use  of  an  average 
value  of  0 then  becomes  physically  unrealistic.  More 
important  is  the  possibility  of  dual  (or  “split")  flow 
due  to  the  change  in  direction  of  the  buoyancy  force. 
Another  factor  is  the  change  of  phase  at  the  solid/ 
liquid  interface,  which  introduces  an  interfacial 
velocity  and  may  distort  the  flow  field,  thereby  affect- 
ing the  heat  transfer  rate.  Tkachev*  was  the  first  to 
observe  the  minimum  Nusselt  number  at  5.5°C  in  his 
experimental  work  with  melting  cylinders.  He  attributed 
this  minimum  to  the  possibility  of  two  opposing  flow 
regimes  within  the  boundary  layer.  Schechter  and 
Isbin*  experimentally  demonstrated  the  existence  of 
dual  flow.  From  their  measurements  of  the  velocity 
profiles  near  a warm  plate  immersed  in  water,  they 
showed  that,  at  certain  temperatures,  the  inner  portion 
of  the  boundary  layer  moved  upwards  while  the  outer 
portion  moved  downwards.  Merks  solved  the  problem 
by  using  an  integral  method  with  a moving  coordinate 
system  and  by  considering  water  as  a high  Prandtl  num- 
ber fluid.  For  the  case  of  ice  spheres,  Merk’s  analysis 
indicated  that  there  was  a minimum  Nusselt  number 
and  a complete  inversion  of  flow  at  a bulk  tempera- 
ture of  4.8°C. 


EXPERIMENTAL  APPARATUS 
General  considerations 

By  the  nature  of  this  experiment,  there  is  a need  for 
the  maintenance  of  a constant  water  temperature  and 
a minimum  heat  loss  in  pumping  water  at  a desired  rate 
through  the  flowmeter  to  the  inlet  of  the  test  assembly. 
PVC  pipe  was  selected  for  use  in  all  piping  because  of 
its  low  rhermal  conductivity.  Clear  Plexiglas  was  used 
for  the  construction  of  the  test  tank  because  of  its 
optica!  clarity  and  easy  fabrication. 

The  other  factor  to  be  considered  was  the  mainte- 
nance of  water  quality  at  constant  composition 
throughout  the  experiments.  The  high  corrosion  re- 
sistance of  aluminum,  stainless  steel,  and  brass  made 
their  use  practical  in  various  load  bearing  frames  and 
components  that  were  in  contact  with  the  water.  A 
small  piece  of  magnesium  was  placed  within  the 
reservoir  as  a sacrificial  metal  to  protect  the  aluminum 
from  cathodic  reaction  with  the  brass  and  stainless 
steel.  Distilled  water  was  used  to  eliminate  the  problem 
of  chlorine  corrosion  found  in  tap  water  and  the  pro- 
duction of  metal  salts  that  could  influence  the  melting 
rate  of  icc.  Designs  using  Teflon  and  nylon  bushing  and 
bearing  surfaces,  as  welt  as  close  sliding  fits  between 


associated  equipment  for  easy  and  accurate  assembly 
of  the  apparatus,  were  specified  to  eliminate  the  need 
for  organic  lubricants  which  could  contaminate  the 
water  and  affect  its  physical  properties.  Apparatus  com- 
ponents conuiningoil,  such  as  gear  reducers  and  flexi- 
ble shafts,  were  separated  from  the  rest  of  the  system. 
This  precluded  the  possibility  of  introducing  detergent 
oils  that  would  combine  with  any  trace  metal  salts  in 
the  water  to  create  detergent  soaps,  which  would  have 
altered  the  viscosity  and  surface  tension  of  the  water 
and  introduced  unknown  factors  into  the  melting 
characteristics  of  the  test  ice  sample. 

Figure  I shows  the  schematic  of  the  experimental 
setup.  It  essentially  consists  of  four  major  parts,  i.e. 
the  water  reservoir,  test  chamber  assembly,  and  pump- 
ing and  control  systems.  The  entire  assembly  was  de- 
signed to  accommodate  a test  sample  of  ice  with  an 
exposed  vertical  face  0.914  m in  height  and  0.30S  m in 
width. 

Water  reservoir 

The  1 .22  x 1 .22  x 0.46  m high  tank  is  constructed 
of  welded  4.76-mm  aluminum  plate.  It  is  supported 
above  the  floor  by  a channel  aluminum  footing  to  re- 
duce the  heat  transfer  between  the  bottom  of  the  tank 
and  the  coldroom  floor.  Both  the  bottom  and  side  of 
the  tank  are  insulated  externally  by  thick  Styro- 

foam sheeting.  The  tank  has  an  operating  capacity  of 
0.473  mJ  and  is  equipped  with  a 7*  horsepower  mixer. 
This  mixer,  through  the  aid  of  the  baffle  system  within 
the  reservoir,  draws  water  past  the  seven  750-W  heaters 
on  the  two  adjoining  sides  of  the  tank,  pulling  it  through 
a passage  in  the  baffle  plates  into  the  major  open  section 
of  the  reservoir  where  it  is  thoroughly  mixed.  A therm- 
istor probe,  placed  within  this  area,  is  connected  to  an 
on-off  controller  for  the  heaters.  This  probe  registers 
temperature  change  to  within  ± 0.05°C. 

Test  chamber  assembly 

The  test  chamber  is  constructed  of  19-mm  clear 
Plexiglas  with  its  external  dimensions  of  0.495  x 0.445 
x 1 .238  m high.  It  contains  the  ice  frame,  ice-frame 
baffles,  internal  frame  belt  assembly,  drain  system  and 
inlet  flow  distribution  box.  Figure  2 shows  the  ice-frame 
and  the  baffle  device  which  has  the  general  dimensions 
of  0.425  m wide  x 0.914  m long  x 0.1 52  m thick.  The 
ice  frame  is  also  constructed  of  19-mm  Plexiglas  with 
8.5-mm  aluminum  plate  mounted  in  back.  The  ice  sam- 
ple has  a surface  area  of  0.914  x 0.305  m exposed  to 
melting.  The  weight  of  the  ice  frame  is  about  19.96  kg 
and  the  ice  sample  weighs  about  42.64  kg.  The  ice-frame 
slips  down  within  recessed  surfaces  on  the  inside  of  the 
test  chamber  and  rests  upon  a baffle  arrangement  at  the 
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1.  Rotameter  6.  Heaters  11.  Internal  frame  belt  assembly 

2.  Mixer  motor  7.  Test  chamber  drain  valves  12.  Belt  drive  assembly 

3.  Flow  control  valves  8.  Control  box  13.  Test  chamber 

4.  Pump  assembly  9.  Ice  face  14.  Thermocouple  tubes 

5.  Reservoir  10.  Overflow  pipe 


Figure  1.  Schematic  of  the  experimental  setup. 


bottom  of  the  chamber,  which  reduces  unwanted 
erosion  on  the  end  of  the  ice  sample.  A similar  arrange- 
ment on  the  top  end  of  the  ice  tramc  serves  the  same 
purpose.  Uniform  flow  is  achieved  by  evenly  distribut- 
ing and  draining  the  water  over  the  entire  cross  section, 
and  simulated  water  flow  on  both  sides  of  the  chamber 
is  accomplished  by  the  use  of  belts  (made  of  plastic 
sheeting)  moving  down  at  the  same  velocity  as  the 
water  flowing  past  the  ice. 

The  flow  distributor,  constructed  of  Plexiglas  with 
a perforated  bottom,  sits  on  top  of  the  test  chamber; 
it  has  0.152-m-high  sides  to  control  splashing  and  is 
fed  by  a four-outlet  header  pipe.  The  water  is  intro- 
duced at  right  angles  to  the  bottom  for  even  water  dis- 
tribution over  the  top  perforated  plate.  Spacers  are 
used  to  vary  the  distance  between  the  perforated 


bottom  and  the  removable  upper  perforated  plate. 

Both  plates  have  12.7-mm-diam  drilled  holes  arranged 
in  a mutually  staggered  fashion  so  that  no  flow  goes 
directly  through  both  plates  without  being  diverted  at 
right  angles  to  reach  the  hole  of  the  bottom  plate.  The 
net  effect  is  to  provide  an  even  watei  flow  throughout 
the  entire  cross  section  of  the  test  chamber.  Similarly 
the  drain  system  is  designed  to  provide  an  even  with- 
drawal of  water  throughout  the  cross  section.  This  was 
accomplished  by  installing  six  14-in.  I.D.  PVC  pipes, 
evenly  spaced,  running  horizontally  out  both  sides  of 
the  test  tank.  A slit,  0.305  m long  and  25.4  mm  wide, 
was  cut  in  each  pipe  in  the  side  facing  the  chamber  bot- 
tom, and  1 2 PVC  gate  valves  were  installed  at  the 
ends  of  the  pipes  for  adjustment  for  equal  flow  rates 
out  of  each  pipe. 

The  internal  frame  belt  assembly  is  situated  within 
the  perimeter  of  the  test  chamber.  The  belts,  made  of 

0. 089-mm-thick  Mylar,  are  set  up  such  that  there  exists 
exactly  a 0.305  x 0.305-m-square  cross  section  for  the 
water  flow.  The  belts  are  driven  downwards  (inside  the 
belt  frame)  at  the  velocity  of  the  water  flow  to  simulate 
an  infinite  medium  of  water  flow  on  both  sides  of  the 
test  chamber.  This  assembly  is  driven  by  a '/u  -horse- 
power variable  speed  d.c.  motor  with  a speed  range  of 
100  to  1 750  rpm.  The  speed  is  reduced  by  a 10:1  ratio 
gear  reducer  with  a single  takeoff  which  in  turn  drives 

a 20:1  ratio  gear  reducer  with  two  takeoffs.  Two  flexi- 
ble shafts  connect  the  final  gear  reducer  to  the  belt 
roller  of  38-mm  diameter.  The  belt  drive  rollers  are 
capable  of  operating  anywhere  within  the  range  of  0.5 
to  8.75  rpm. 

The  entire  test  chamber  and  belt  drive  assembly  is 
supported  by  the  external  frame  centrally  situated  in 
the  reservoir.  This  frame  is  mounted  on  threaded  brass 
studs  on  the  reservoir  bottom  and  can  be  adjusted  such 
that  the  entire  frame  and  assembly  are  vertical.  The 
test  chamber  is  hinged  on  the  external  frame  so  that  it 
may  be  tilted  back  for  easy  removal  of  the  belt  drive 
assembly,  ice  frame  and  samples. 

Pumping  system 

The  reservoir  water  is  pumped  by  a rotary  gear  pump 
operating  at  its  maximum  rate  of  4.52x1 0J  kg/h.  The 
pump  assembly  is  connected  to  the  reservoir  by  a flexi- 
ble reinforced  rubber  hose.  This  reduces  the  vibration 
that  would  have  been  transmitted  to  the  tank  if  a solid 
pipe  had  been  used,  eliminating  considerable  noise  and 
associated  fatigue  of  welds.  The  desired  amount  of  water 
discharged  to  the  test  chamber  is  regulated  by  a bypass 
and  a main  line  flow  control  system.  It  consists  of  1-in. 
PVC  pipe  with  the  main  line  flow  controlled  by  a 1-in.- 

1. D.  brass  gate  valve  leading  to  the  Rotameter  and  a 
bypass  line  upstream  of  the  valve  returning  the  excess 
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water  to  the  reservoir.  The  Rotameter  is  rated  from 
1.388X10*4  to  1.388x10* 3 m3/s. 

Electrical  control  system 

All  electrical  controls  for  the  experiment  are  arranged 
together  within  a control  box  as  an  integral  part  of  the 
experimental  setup.  Both  1 20-  and  208-V  power  is 
supplied  through  individual  disconnect  switches.  A 
variable  speed  control  unit  is  used  to  regulate  the  ‘/n- 
horsepower  d.c.  motor  for  the  belt  drive.  A thermistor 
on-off  control  connected  to  the  probe  in  the  reservoir 
tank  operates  three  relays  for  the  seven  heaters  in  the 
reservoir.  The  switch  panel  provides  independent 
manual  control  for  each  of  the  seven  heaters,  as  well  as 
the  mixer  motor,  rotary  gear  pump,  and  temperature 
control  unit. 


EXPERIMENTAL  PROCEDURE 

The  ice  samples  are  made  by  first  preparing  the  ice 
frame  for  freezing  the  distilled  water.  During  this 
process,  the  end  plates  of  12.7-mm-thick  Plexiglas  are 
attached  to  both  ends  of  the  ice  frame.  The  plates  are 
held  in  place  by  screws  with  all  mating  surfaces  between 
the  ice  frame  and  end  plates  being  temporarily  sealed 
with  silicone  elastomer.  The  ice  anchor,  made  of  two 
PVC  tubings  leading  to  threaded  holes  in  a block  of 
Plexiglas,  then  is  put  into  a position  within  the  ice  frame 
with  the  ends  of  the  two  tubings  plugged  with  rubber 
stoppers  and  also  sealed  with  silicone  elastomer  to  pre- 
vent water  from  entering  (see  Fig.  2).  The  anchor  is 
kept  from  moving  out  of  place  by  lead  weights,  and  the 
plugged  ends  of  the  tubings  rest  against  the  end  plate. 
The  ice  frame  is  then  placed  and  leveled  on  the  alumi- 
num cold  plates;  this  enables  the  ice  surface  to  be  ievel 
with  the  frame’s  front  surface.  Ethylene  glycol  at 
-20°C  circulates  through  the  0.28-m-diam  coldplates 
from  a constant  temperature  bath  cooled  by  trichloro- 
ethylene at  -73.0°C.  In  filling  the  ice  frame  with  dis- 
tilled water,  the  expansion  factor  in  forming  ice  from 
water  is  also  taken  into  account.  During  the  entire 
freezing  period,  the  ice  frame  is  covered  with  an  in- 
verted, tightly  fitting  38-mm-thick  Styrofoam  box 
(without  top)  and  the  coldroom  temperature  is  kept  at 
1.67°C.  This  allows  the  ice  to  form  in  an  upward 
direction  only  and  eliminates  the  ice  expansion  problems 
associated  with  the  phase  transition.  Once  the  ice  sam- 
ple is  made  (which  usually  takes  four  days)  the  surface 
is  smoothed  and  evened  with  a warm  copper  billet.  The 
excess  ice  is  removed,  making  the  ice  surface  flush  with 
the  ice  frame.  The  ice  frame  is  then  placed  in  a cold- 
room  maintained  at  -1 ,0°C  and  remains  there  for  at 
least  36  hours  with  a fan  blowing  on  the  backplate  to 


ensure  that  the  frame  and  the  ice  are  conditioned  to  a 
uniform  temperature  of  -1.0 °C. 

Prior  to  the  experiment,  the  end  plates  and  plugs  for 
the  ice  anchor  are  removed  and  the  silicone  elastomer 
cleared  from  the  ice  frame.  Then  two  stainless  steel  rods 
are  passed  down  through  the  tubes  of  the  ice  anchor  and 
threaded  into  the  Plexiglas  block  to  create  a handle  for 
placing  the  ice  frame  in  and  removing  it  from  the  test 
reservoir. 

In  the  meantime,  the  apparatus  is  being  prepared  for 
the  experiment.  Before  the  test,  the  12  drain  valves  are 
calibrated  for  a specific  flow  velocity  by  volumetric 
flow  timing  with  each  one  adjusted  to  equal  flow  rates 
at  the  test  temperature.  Then  the  coldroom  tempera- 
ture is  raised  to  20°C  to  dry  out  the  water  remaining  in 
the  valve  seats  and  other  parts  of  the  piping  system  to 
prevent  possible  freezeup  when  the  room  temperature  is 
brought  down  to  -1 .0°C  during  the  experiment.  To 
assure  a plug  flow  in  the  test  chamber,  the  variable  speed 
drive  is  adjusted  so  that  the  timed  movement  of  the  belts 
matches  the  calculated  water  velocity  in  the  chamber. 

The  temperature  for  the  test  is  set  by  bringing  the 
reservoir  to  the  desired  temperature  and  maintaining  it 
at  that  level  with  the  on-off  thermistor  control. 

Prior  to  the  experiment,  the  coldroom  temperature  is 
adjusted  to  that  of  the  conditioned  ice  sample  (-1 ,0°C). 
With  the  reservoir  mixer  running,  the  heater  control  sys- 
tem is  turned  on  to  establish  the  desired  steady  state 
temperature  in  the  reservoir.  Once  the  room  is  at 
-1 .0°C,  the  ice  frame  is  brought  to  the  room  and 
weighed,  the  test  tank  tilted  back,  and  the  ice  frame 
placed  into  the  chamber.  The  handle  attached  to  the 
anchor  is  removed  and  replaced  with  two  rubber  stoppers, 
and  the  internal  frame  belt  assembly  is  slipped  into  the 
chamber.  The  inlet  flow  distribution  box  is  attached  to 
the  side  wall  of  the  test  chamber  and  two  thermistors 
are  then  inserted  through  the  perforated  plate  and  left 
hanging  more  or  less  one-third  and  two-thirds  down  the 
ice  frame  to  measure  the  variation  of  temperature  in  the 
water  during  the  experiment  (the  original  design  with 
thermocouples  was  not  used).  A 1 ,6-mm  Plexiglas  sheet 
is  inserted  between  the  belt  assembly  and  ice  frame  (to 
prevent  melting  during  rapid  filling  of  the  chamber)  and 
the  top  ice  frame  baffle  is  installed.  The  test  chamber  is 
then  erected  into  vertical  position  and  the  inlet  header 
pipe  and  flexible  drive  shafts  are  connected. 

An  experiment  is  initiated  by  starting  the  pump  and 
closing  the  bypass  valve.  The  belt  drive  is  turned  on 
while  the  test  chamber  is  being  filled.  As  soon  as  the 
water  reaches  the  overflow  level,  the  Plexiglas  sheet  is  re- 
moved and  the  experiment  started.  The  flow  rate  is 
adjusted  by  use  of  the  mainline  and  bypass  valves  until 
there  is  a small  trickle  coming  out  the  overflow  pipe. 

The  desired  flow  rate  in  the  test  chamber  is  regulated  by 
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Figure  3.  Melting  rate  as  a function  of  bulk  water  temperature. 


setting  the  Rotameter  reading.  The  thermistor  readings  ice  melted,  the  variation  of  the  extent  of  ice  melting  is 
are  recorded  periodically  to  see  if  any  variations  in  measured  at  various  points  on  the  surface, 

temperature  existed  during  the  test.  The  duration  of 
the  experiment  is  determined  primarily  by  the  water 

temperature.  For  temperatures  as  low  as  1 .1°C,  experi-  EXPERIMENTAL  RESULTS 
ments  last  as  long  as  six  hours  and  still  do  not  have  an 

excessive  amount  of  melting.  For  higher  water  tempera-  In  total,  46  experiments  were  conducted  covering 

tures,  the  duration  is  usually  about  two  or  three  hours.  test  chamber  water  velocities  v of  1 .7  to  9.8  mm/s  and 

The  experiment  is  terminated  by  turning  off  the  gear  bulk  water  temperatures  7.  of  1 .11°  to  7.50°C.  The 

pump  and  unplugging  the  two  auxiliary  drains  located  ice  sample  prior  to  the  experiment  had  a smooth  surface 

at  the  bottom  of  the  test  chamber.  Then  the  inlet  flush  with  the  ice-frame,  but  after  the  experiment,  the 

water  pipe  is  disconnected,  and  the  thermistor  leads  ice  surface  was  found  to  be  relatively  flat  with  some 

removed  from  the  chamber.  The  test  chamber  is  tilted  extra  melting  along  the  side  near  the  top.  This  is 
back  as  soon  as  the  water  level  reaches  the  bottom  of  attributed  to  the  turbulence  created  by  the  inside  edges 

the  chamber  wd  the  inlet  flow  distributor,  top  ice  of  the  ice-frame  as  ice  melts  away.  In  most  cases,  with 

frame  baffle  and  belt  drive  assembly  are  removed.  the  exception  of  the  highest  water  velocity  used  in  the 

Since  the  ice  and  the  room  are  maintained  at  -1 ,0°C,  experiment,  the  ice  surface  at  the  termination  of  the 

a thin  layer  of  ice  forms  between  the  backplate  of  the  experiment  remains  flat  indicating  even  melting.  For 

ice  frame  and  the  test  chamber  during  the  test.  To  the  case  of  9.8-mm/s  velocity  and  the  highest  tempera- 

ease  removal  of  the  ice  frame  from  the  chamber,  usually  ture  (i.e.  7.50°C)  there  is  a greater  melting  near  the 

a small  amount  of  hot  water  is  poured  directly  on  this  water  entrance  to  the  test  chamber.  Experiments  1 2-1 3 

interface.  As  soon  as  the  ice  frame  is  removed  from  the  and  14-15  were  run  to  determine  the  reproducibility  of 

test  chamber  and  weighed  to  determine  the  amount  of  data  under  identical  conditions.  It  can  be  noted  that 
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Table  I.  Experimental  parameters  and  computed  results. 


Experiment 

number 

Bui * water 
temperature 

Tmrc) 

Experiment 

duration 

»(h) 

Water 

velocity 

v(mm/$) 

Melting 

rate 

mfkg/h) 

Husselt 

number 

NU|_ 

Prandtl 

number 

fr 

Reynolds 

number 

r 

1 

2.2 

6.00 

5.1 

1.60 

396 

13.0 

3.03 

2 

1.94 

3.00 

5.1 

1.21 

342 

13.1 

2.80 

3 

2.22 

6.00 

1.7 

1.07 

264 

13.0 

0.87 

4 

2.22 

6.00 

1.7 

1.20 

296 

13.0 

0.94 

5 

4.72 

4.00 

1.7 

2.09 

241 

12.4 

1.02 

6* 

4.72 

4.00 

1.7 

2.25 

259 

12.4 

1.02 

7 

1.94 

5.00 

9.8 

1.49 

419 

13.1 

6.00 

8 

1.94 

5.00 

9.8 

1.60 

452 

13.1 

6.03 

9* 

1.94 

5.00 

9.8 

1.69 

478 

13.1 

5.65 

10 

1.83 

5.00 

9.8 

1.47 

441 

13.1 

5.53 

11* 

1.83 

5.00 

9.8 

1.53 

462 

13.2 

5.46 

12 

7.50 

1.25 

9.8 

6.11 

443 

11.8 

7.07 

13 

7.50 

1.25 

9.8 

6.10 

442 

11.8 

7.07 

14 

4.72 

2.50 

9.8 

3.76 

437 

12.4 

6.50 

IS 

4.72 

2.50 

9.8 

3.67 

424 

12.4 

6.50 

16 

2.77 

4.00 

9.8 

2.07 

408 

12.6 

6.10 

17 

3.61 

3.00 

9.8 

2.75 

416 

12.7 

6.64 

18 

5.83 

2.00 

9.8 

4.65 

433 

12.2 

7.32 

19 

5.83 

2.50 

5.1 

3.68 

344 

12.2 

3.40 

20 

7.20 

2.00 

5.1 

4.64 

350 

11.9 

3.54 

21 

4.44 

3.00 

5.1 

2.95 

362 

11.9 

3.25 

22 

3.33 

4.50 

5.1 

2.10 

343 

12.3 

3.15 

23 

1.11 

7.00 

5.1 

0.76 

377 

13.3 

2.91 

24 

5.83 

2.50 

1.7 

2.92 

290 

12.2 

1.13 

25 

7.22 

2.00 

1.7 

3.88 

308 

11.9 

1.18 

26 

4.44 

3.00 

1.7 

2.52 

329 

12.5 

1.12 

27 

5.56 

1.75 

1.7 

2.76 

268 

12.3 

1.12 

28 

5.80 

2.75 

1.7 

2.71 

252 

12.2 

1.17 

29 

3.47 

5.00 

1.7 

1.79 

280 

12.7 

1.10 

30 

2.42 

5.00 

1.7 

1.02 

227 

13.0 

1.05 

31 

2.25 

5.00 

1.7 

0.92 

222 

13.0 

1.05 

32 

3.84 

3.00 

1.7 

1.83 

258 

12.6 

1.10 

33 

4.41 

3.00 

1.7 

2.33 

285 

12.5 

1.13 

34 

4.95 

3.00 

1.7 

2.43 

265 

12.4 

1.15 

35 

6.18 

2.50 

1.7 

3.27 

284 

12.1 

1.19 

36 

2.25 

4.75 

1.7 

1.14 

277 

13.0 

1.05 

37 

154 

5.00 

1.7 

0.73 

261 

13.2 

1.28 

38 

3.30 

2.50 

6.1 

2.35 

370 

12.8 

4.68 

39 

5.80 

2.17 

6.1 

4.08 

379 

12.2 

3.74 

40 

5.75 

3.00 

3.4 

3.66 

343 

12.2 

2.08 

41 

7.04 

2.58 

3.4 

4.49 

341 

11.9 

2.12 

42 

3.64 

4.42 

3.4 

2.17 

324 

12.7 

1.92 

43 

2.40 

5.00 

3.4 

1.47 

334 

13.0 

1.84 

44 

4.43 

2.50 

5.1 

2.89 

352 

12.5 

2.95 

45 

4.43 

2.50 

5.1 

3.09 

376 

12.5 

2.95 

46 

4.43 

2.50 

5.1 

3.25 

369 

12.5 

2.95 

* Experiments  withe  belts  operating. 


the  apparatus  was  able  to  provide  consistent  data.  A 
few  experimental  runs  were  also  conducted  without 
the  belts  moving  to  assess  the  effect  of  the  buildup  of 
side  boundaries  on  the  melting  rate.  To  simulate  the 
prototype  operation  conditions  more  closely,  all  other 
experiments  were  conducted  with  the  belts  moving. 
Table  I shows  the  experimental  conditions  and  the 


calculated  results.  Figure  3 shows  the  variation  of 
melting  rate  M wi  *h  Tm,  using  v as  a parameter.  All 
the  curves  are  extrapolated  to  point  zero  for  no  melting 
at  0°C.  The  data  show  a great  deal  of  scattering, 
demonstrating  the  possibility  that  varying  extents  of 
turbulence  were  introduced  at  the  entrance  of  the  test 
chamber.  As  mentioned  earlier,  the  water  must  be 
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maintained  at  an  overflow  level  during  the  entire  test 
period.  Otherwise,  there  will  be  a gap  between  the 
inlet  water  distributor  and  the  water  level  in  the  test 
chamber,  and  the  water  will  be  running  down  to  the 
chamber,  thus  creating  turbulence  and  promoting  a 
higher  rate  of  heat  transfer.  Also  the  low  sensitivity 
of  the  weighing  scales  has  contributed  to  the  incon- 
sistency of  the  data.  Regardless  of  the  scattering  of 
the  data,  however,  the  trends  clearly  show  that  the 
ice  melting  rate  has  a strong  dependence  on  bulk  water 
temperature  and  water  velocity. 

The  experimental  results  are  expressed  in  terms  of 
an  average  heat  transfer  coefficient  h defined  as 


where  qi  - sensible  heat  required  to  raise  the  ice  from 
initial  temperature  (<  0°C)  to  the  melt- 
ing temperature  Tm  (=  0°C) 

<7f  = heat  of  fusion  associated  with  melting 
A = area  of  the  melting  surface 
AT  = temperature  difference  between  Tm  and 
Tm  (since  Tm  =0°C,  AT=  Tm) 

0 = the  time  duration  of  the  experiment 


A Reynolds  number  (ReL  = vpL/p)  based  on  ice  plate 
height  is  used  to  define  the  flow  characteristics,  where 
p and  n are  the  water  density  and  viscosity,  respectively. 
The  physical  properties  are  evaluated  at  the  arithmetic 
mean  of  bulk  water  temperature  Tm  and  melting  tem- 
perature Tm  or  TJ2.  For  each  experiment,  the  weight 
of  ice  melted  W,  as  well  as  Tm,  6 and  the  flow  rate  are 


recorded.  The  amount  of  ice  melted  W is  used  to  com- 
pute q,  = - r,),andqf  = WLf  where  c„  T, 

and  Lf  are  the  specific  heat,  initial  temperature,  and 
latent  heat  of  fusion  of  the  ice,  respectively.  Figure  4 
shows  a correlation  between  NuL/Pr  and  ReL  in  which 
Nul  and  ReL  are  evaluated  based  on  height  of  the  plate. 
The  relation  can  be  fairly  represented  by 

^ = 3.27 S (ReL)°-270  (2) 

with  a correlation  coefficient  of  0.843.  The  relatively 
low  value  of  the  correlation  coefficient  is  essentially 
due  to  the  difficulty  in  ascertaining  the  accuracy  in 
the  temperature  measurement,  especially  when  the  bulk 
water  temperature  is  near  the  melting  temperature.  For 
AT  = 1°C,  a 0.1°C  error  introduces  a 10%  variation  in 
the  heat  transfer  coefficient,  assuming  that  the  ice 
weighing  is  reproducible.  A few  experimental  data 
without  the  belts  moving  are  also  indicated  in  Figure  4. 
It  can  be  seen  that,  although  there  seems  to  be  some 
effect  of  the  moving  belts  on  the  ice  melting  rate,  the 
effect  of  the  two  side  plates  (belts)  does  not  reduce  the 
heat  transfer  rate  noticeably,  and  in  fact  the  data 
are  all  within  the  experimental  range  from  one  to 
another.  A partial  explanation  is  that  there  would  be 
minimal  frictional  resistance  to  flow  with  a very  thin 
boundary  layer  buildup  because  the  belts  are  made  of 
smooth  Mylar  sheets.  Another  explanation  is  that  the 
ice  plate  is  wide  enough  to  consider  the  boundary  effect 
insignificant. 
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DISCUSSION  AND  COMPARISON  OF  RESULTS 


Therefore,  the  Nuswlt  number  in  caw  of  melting  is 


r 

■ i 
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As  indicated  in  the  Introduction,  there  have  been 
no  experimental  forced  flow  studies  of  a vertical 
melting  plate.  However,  Yen  and  Tien11  solved,  in  an 
analytical  study,  the  laminar  heat  transfer  problem  over 
a horizontal  melting  plate  with  a linear  boundary  layer 
velocity  distribution,  i.e.  vx  - By  where  vx  is  the 
velocity  along  the  plate.  The  equation  of  energy  was 
solved  analytically  and  the  Nuswlt  number  can  be 
expreswd  by 


where  L is  plate  length.  The  value  of  an  is  evaluated 
from  the  integral 

an  = f exp(-A3  + 0A/<7n)<A  (4) 

0 

where  0 = c(7\.  - Tm)ILf  and  can  be  regarded  as  a 
melting  parameter.  The  value  of  an  is  found  to  be  a 
function  of  0.  By  approximating  the  value  of  B with 
rwln , where  rw  is  ice-water  interface  shear  stress  and 
H is  the  dynamic  viscosity  of  water,  the  ratio  of  NuL 
(with  melting)  to  Nu0  (without  melting)  is  found  to 
be 

NuL  _ £o  AwY/3 
Nuo  on  \^0  J 

where 

m 

aQ  = / exp  (-\3)<A. 

0 

Since  the  velocity  profile  is  assumed  to  be  com- 
parable to  the  temperature  profile  and  r0  is  pro- 
portional to  the  velocity  gradient,  eq  S can  be  reduced 
to 


(5) 


(6) 


For  the  caw  of  flow  over  a flat  plate  with  no  pressure 
gradient  and  no  mass  transfer  across  the  plate,  the 
Nuswlt  number  is  given  by  Schlichting7  as 

Nu0  = 0.664  (Pr)1/3  (ReL)3/2  . (8) 


NuL  = 0.664  <Pr)1/3  (Re|.),/J  • W 

Figure  5 shows  a comparison  of  eq  8 and  9 and  the 
experimental  data  of  the  prewnt  investigation.  Experi- 
mental conditions  are  uwd  to  evaluate  the  value  of  the 
parameter  ( o0/on )4^3  in  eq  9,  which  depends  on 
c TJLf.  Since  c TJLi  varies  in  the  experiments  be- 
tween 0.014  and  0.094,  the  effect  of  (a^/an)*^  on  the 
values  of  NuL  is  rather  small.  However,  for  c Tm/Lf  = 

1 .0  the  reduction  in  heat  transfer  can  be  quite  large, 
approximately  amounting  to  42%.  As  can  be  noted 
from  the  figure,  the  prewnt  experimental  results  arc 
about  four  to  five  times  greater  than  thow  predicted 
from  a horizontal  flat  melting  plate.  The  results  can  be 
fairly  reprewnted  by 

= 16.5(ReL)°-279  (10) 

/is  tl/2  L 


with  a moderate  correlation  coefficient  of  0.89S. 

A few  data  points  from  Stubstad  and  Quinn*  are 
shown  in  Figure  4.  The  disagreement  with  the  present 
data  is  apparent.  In  their  experiment,  a cylindrical 
column  of  ice  3 m in  height  and  of  1 .83  m initial  diameter 
was  uwd.  For  a given  mass  flow  rate,  the  annular  space 
between  the  tank  wall  and  ice  surface  increases  as  melt- 
ing progresws;  therefore,  for  a given  experiment  the 
value  of  ReL  is  reduced.  In  their  computation  of  an 
average  heat  transfer  coefficient,  the  assumption  was 
made  that  there  was  no  change  in  height  during  the 
experiment.  (In  reality,  this  is  not  true.)  The  value  of 
h is  computed  from 


where  Tm  is  average  bulk  water  temperature.  The  value 
of  incremental  change  of  ice  radius  R with  time,  ob- 
tained by  considering  the  volume  and  mass  change  of 
ice,  can  be  finally  expreswd  as 

A R = /?«?,)  -V*2(0, ) - (A0/p,irL)(<ttf,/tf»)  (12) 

where  M,  is  the  initial  ice  mass  in  the  tank.  From  this 
relationship  it  is  possible  to  determine  the  average 
change  of  the  radius  of  the  ice  cylinder  during  any  time 
interval  from  0,  to  02.  The  diwrepancy  between  thew 
data  and  the  prewnt  data  can  be  attributed  largely  to 
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Figure  5.  Comparison  of  present  data  with  computed  values  from 
horizontal  plate. 


the  assumption  of  constant  height  and  the  method  of 
obtaining  the  variation  of  R with  the  water  level  change 
associated  with  the  phase  change.  The  higher  values 
of  Nul  are  also  due  to  the  high  degree  of  turbulence 
of  the  entering  water  spraying  directly  down  from  the 
nozzle  into  the  annular  gap. 

As  mentioned  in  the  Introduction,  quite  a few  in- 
vestigators have  worked  on  the  free  convective  melting 
heat  transfer  problem.  Vanier  and  Tien10  studied  the 
problem  both  analytically  and  experimentally.  By 
studying  a combination  of  plate  and  bulk  water  tem- 
peratures, they  were  able  to  identify  the  various  flow 
regions  created  by  the  unique  property  of  water  which 
exhibits  a maximum  density  at  about  4°C.  With 
boundary  layer  approximations,  including  the  effect 
of  density  inversion,  and  the  use  of  similarity  trans- 
formation, they  solved  the  equations  of  motion  and 
energy  numerically.  The  Nusselt  number  was  found  and 
expressed  as 

Nul  =2^I(Gr),/4|_^(0)j  (13) 

where  the  Grashof  number  is  always  positive  and  is 
defined  as 


Gr 


_ g/-3 |g,  (r„  - rj 


where  Tp  is  plate  temperature,  and  v the  kinematic 
viscosity.  0.  is  defined  by 


= 01  +202?-  ♦ 303 

1+0,  Tm+02Tl+fi3Tl 


OS) 


in  which  01 , 02  and  03  are  the  coefficients  in  the  cubic 
expression  of  the  density-temperature  relationship  of 
water  and  H ' (0)  is  the  dimensionless  temperature  gra- 
dient at  the  plate.  For  the  case  of  7"p  = 0°C  (equiva- 
lent to  the  case  of  melting  Tm  = 0°C),  Vanier  and 
Tien’s  numerical  results  (after  scaling  to  the  present 
experimental  plate  height)  are  shown  in  Figure  6. 

The  most  recent  experimental  study  of  natural 
convection  flow  over  a vertical  ice  slab  immersed  in 
cold  water  has  been  reported  by  Bendetl  and  Gebhart.1 
In  this  study,  the  effect  of  the  bulk  water  temperature, 
especially  in  the  vicinity  of  the  density  extremum,  on 
heat  transfer  was  determined.  Also,  the  effect  of  the 
density  maximum  on  flow  direction  and  bulk  water 
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Figure  6.  Comparison  of  present  data  with  those  from  free  con- 
vection studies. 


temperature  where  up  flow  changes  to  downflow  was  The  greater  scattering  of  data,  especially  at  the  lower 

determined.  They  compared  their  experimental  results  extreme  of  the  temperature  used  in  this  experiment, 

with  the  analytical  work  of  Gebhart  and  Mollendorf3  can  be  partially  attributed  to  the  inaccuracy  in  the 

and  found  a mean  magnitude  of  deviation  of  5.6%.  temperature  measurement.  A small  error  in  the 

Their  experimental  results  are  also  shown  in  Figure  6 absolute  temperature  measurement  can  cause  a great 

along  with  the  present  results.  It  can  be  seen  that  the  percentage  of  deviation  and  thus  create  a much  greater 

work  of  Bendell  and  Gebhart  agrees  with  the  analytical  magnitude  of  scattering.  However,  the  effect  of  flow 

work  of  Vanier  and  Tien10  fairly  well.  It  also  can  be  on  the  melting  heat  transfer  can  easily  be  visualized  in 

noted  that,  for  the  case  of  forced  motion,  there  is  a Figures  3, 4 and  5.  For  the  temperature  range  covered 

significant  increase  in  heat  transfer  rate  over  free  con-  in  this  study  it  seems  that  the  change  in  the  direction 

vection  even  with  a minimal  velocity  of  1 .7  mm/s.  of  flow  due  to  buoyancy  forces  has  exhibited  no 
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^ = 3.275  (ReL)*™> 

with  a moderate  correlation  coefficient  of  0.843.  Al- 
though the  experimental  runs  were  only  conducted  for 
one  specific  plate  height,  this  equation  will  probably 
predict  the  actual  heat  transfer  characteristics  of  the 
prototype. 

3.  From  this  study  it  can  be  concluded  that,  in  the 
case  of  forced  convective  melting  heat  transfer,  the 
effect  of  buoyancy  forces  caused  by  the  density  maxi- 
mum at  4°C  on  heat  transfer  appeared  insignificant. 
This  is  true  even  at  very  low  flow  velocity,  as  the  heat 
transfer  is  drastically  reduced  in  the  vicinity  of  5.6°C 
in  the  case  of  free  convective  melting  heat  transfer. 
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significant  effect  on  melting  rate,  as  is  indicated 
numerically  by  Vanier  and  Tien,10  and  experimentally 
by  Bendell  and  Gebhart.' 


CONCLUSIONS 

1 . An  experimental  device  to  simulate  the  rather 
large  dimensions  of  a prototype  heat  sink  has  been 
developed.  This  was  accomplished  by  installing  belts 
on  two  opposite  sides  of  a square  column.  The  belts 
were  designed  and  adjusted  to  move  at  the  same 
direction  and  velocity  as  the  water  flowing  through 
the  test  chamber.  In  doing  this,  the  laboratory  setup 
could  be  considered  as  a segmental  ring  cut  out  from 
the  proposed  prototype. 

2.  For  bulk  water  temperatures  ranging  from  1.11 
to  7.5°C  and  water  velocities  from  1.7  to  9.8  mm/s, 
the  present  results  could  be  fairly  represented  by 
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In  accordance  with  letter  from  DAEN-RDC,  DAEN-ASI  dated 
22  July  1977,  Subject:  Facsimile  Catalog  Cards  for 
Laboratory  Technical  Publications,  a facsimile  catalog 
card  in  Library  of  Congress  MARC  format  is  reproduced 
below. 
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